Evidence for Dirac Nodes from Quantum Oscillations in SrFe 2 As 2 
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We present a detailed study of quantum oscillations in the antiferromagnetically ordered pnictide 
compound SrFe2As2 as the angle between the applied magnetic field and crystalline axes is varied. 
Our measurements were performed on high quality single crystals in a superconducting magnet, 
and in pulsed magnetic fields up to 60 T, allowing us to observe orbits from several small Fermi 
surface pockets. We extract the cyclotron effective mass m* and frequency F for these orbits and 
track their values as the field is rotated away from the c-axis. While a constant ratio of m* /F is 
observed for one orbit as expected for a parabolic band, a clear deviation is observed for another. 
We conclude that this deviation points to an orbit derived from a band with Dirac dispersion near 
the Fermi level. 

PACS numbers: 71.18.+y, 75.30.Fv, 71.20.Ps 



Electrons with Dirac-like dispersion have been invoked 
to explain many fascinating phenomena in materials. 
Arising in systems where the low-energy states are gov- 
erned by the relativistic Dirac equation, they are charac- 
terized by an energy which depends linearly on wavevec- 
tor, and possess high velocities and a zero effective mass. 
Dirac electrons determine the remarkable electronic prop- 
erties of graphene [T] , and in topological insulators they 
may give rise to novel spin textured states through strong 
spin-orbit coupling [2] . Recent experimental work has es- 
tablished a compelling body of evidence for the existence 
of such electrons - in graphene H] and graphite [S] for 
instance, and in topological insulators such as Bi2Se3 [6]. 

In this letter we report evidence from quantum oscilla- 
tion experiments for the existence of Dirac fermions aris- 
ing in SrFe2As2, an antiferromagnetically ordered metal 
which becomes superconducting by the application of hy- 
drostatic pressure [3 Hj or by chemical substitution [9] . 
SrFe2 As2 possesses a Fermi surface of several small pock- 
ets, formed due to the reconstruction of a larger param- 
agnetic Fermi surface by magnetic ordering in the FeAs 
planes along the Q=(ir, n,0) direction [10]. Our study 
tracks the effective mass m* of cyclotron orbits [TT] with 
frequency F that arise from these pockets, as the angle 
between the applied magnetic field and crystalline axes 
is varied. When the orbits enter the polar regions of one 
Fermi surface, we observe a decrease in the ratio m* / F 
from the value expected for a parabolic band. 

We successfully model this decrease in terms of a Dirac 
dispersion at the Fermi level, arising in an effective two- 
band model parameterized by only the interlayer hopping 
t and chemical potential p. This approach was suggested 
previously |12j , but until now complete angular data was 
not available to confirm it. The broad agreement be- 
tween our measurements and this simple model estab- 
lishes the existence of massless excitations and constrains 
their properties. 

Quantum oscillatory studies offer an extremely sensi- 
tive bulk probe of electronic structure. Unlike transport 



measurements, which take into account a weighted aver- 
age of contributions from various bands, quantum oscil- 
lation studies are able to selectively access information 
about each band individually. The presence of such par- 
ticles in our measurements strongly suggests a topologi- 
cally protected nodal spin density wave state [13] , which 
occurs due to a crossing of hole and electron bands near 
the Fermi level resulting from the magnetic reconstruc- 
tion of the Fermi surface. Intriguingly, Dirac fermions 
arising in this manner are proposed to be fundamentally 
different than those in graphene. In the pnictide case, 
the presence of non-Dirac electrons leads to two Dirac 
cones with the same chirality, with distinct implications 
for transport properties [13] . 

Our samples were grown using a self-flux technique 
[T5] , yielding high quality platelet-like single crystals with 
dimensions of approximately 1.5 x 1.5 x 0.05 mm 3 that 
were free of flux inclusions. To improve sample qual- 
ity, single crystals were annealed at 800° C for 24 hours 
under argon gas, with subsequent annealing for 2 hours 
at 300° C in some instances to remove a superconduct- 
ing signal near 20 K [16] without otherwise altering the 
Pt curve. Annealed single crystals were of high qual- 
ity, with P3oqk/P4:K values of up to 17 being reached, 
Over 20 annealed samples were then screened for best 
quantum oscillations. Quantum oscillations at low fields 
down to 12T (Fig. la) were measured for a crystal with 

P300k/P4K ~ 9. 

Oscillations were studied using two techniques, via re- 
sistivity measurements conducted on a superconducting 
cryomagnetic system in Cambridge, with fields as high 
as 16 T, and with a contactless tunnel diode oscillator 
(TDO) technique [IS] using the 60 T short pulse mag- 
net at Los Alamos National Laboratory at temperatures 
from 0.5 K to 5 K. In both systems we are able to vary the 
orientation of the crystal by use of a mechanical rotator. 

Figure 1 shows examples of our data, with 9 defined 
as the angle between the applied magnetic field and the 
crystalline c-axis. The top panel shows results of a four- 
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FIG. 1: (color online) Quantum oscillations in SrFe2As2, plot- 
ted against inverse field. Panel (a) shows resistivity data 
taken in a superconducting magnet, while panels b) and c) 
show the oscillatory component of the resonance frequency of 
a TDO circuit in a pulsed magnetic field. Data is shown for 
different orientations of the field to the crystalline c axis, de- 
fined by 9. The solid lines are a fit to a two-frequency model, 
described in the text. The insets show the Fourier spectra, 
dominated by two peaks which we identify as F a and Fp. 



wire resistivity measurement in a 16 T superconducting 
magnet, while the bottom panels show the resonance fre- 
quency of a TDO circuit which was coupled inductively 
to a sample subjected to a pulsed magnetic field. In both 
cases oscillations arise due to changes in the magnetore- 
sistance of the crystal, described by the Shubnikov-de 
Haas (SdH) effect. A fourth order polynomial has been 
used to subtract the field-dependent background, which 
arises due to the magnetoresistance of the sample. Our 
data is plotted in units of 1/B, indicating an oscillatory 
component of the signal that is periodic in inverse field. 

At the lowest angles, we observe oscillations down to B 
~ 12 Tesla, whose amplitude is exponentially dependent 
on field. The data is well described by a simple fit to a 
model with two frequencies, F a and Fp, with amplitudes 
A and phases S. The total oscillation amplitude is then 
A = A a sin(27rF a /B + 5 a ) + Apsin(2nFp/B + Sp), and 
the results of this fit are shown by the solid lines through 
the data. A Fourier transform of each data set is shown 
in the insets of Fig. 1, showing a spectrum dominated by 
two frequencies which are evident at all angles. 

We identify these frequencies as corresponding to or- 
bits about small pockets of the reconstructed Fermi sur- 
face. In keeping with convention for the pnictide '122' 
compounds [TUl US IE] we label these as (3 and a for the 
lower and higher frequencies respectively, with Fa = 160 
T and F a = 400 T for 6>=10°. At higher angles the over- 
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FIG. 2: (color online) Amplitude of oscillations as a function 
of temperature, defined by the height of the peaks in the 
Fourier spectrum. The amplitudes are normalized to the low 
temperature values. The lines through the data represent fits 
to the Lifshitz-Kosevich damping factor used to extract m*. 



all amplitude of the oscillations is reduced, however we 
are still able to resolve the a and /3 frequencies. For the 
highest angles, a third frequency, F=200 T at 6 = 74° 
becomes resolvable, which we ascribe to the 7 orbit. 

The magnetic field dependence of the amplitude for the 
0=10° data yields Dingle temperatures of To a — 0.7 K 
and Tjjp = 2.7 K. From these we deduce mean free paths 
of 2700 A and 800 A respectively, using Fermi surface 
averaged velocities of v a — \ / 2heF a /m k — 5 x 10 4 ms^ 1 
and V/3 — 6 x 10 4 ms _1 . These long mean free paths are 
comparable to those reported in BaFe2As2 [TO] . 

In Fig. 2 we show how the amplitude of the oscillations 
depends on temperature. We plot the height of the peak 
in the Fourier spectrum as a function of temperature, 
normalized to the low temperature value. This data is 
then fit with the Lifshitz-Kosevich (LK) thermal damp- 
ing factor, R T = X/sinh(X) where X=14.69 m*T/B 
|11) . For the high field data set, we worked with Fourier 
windows between 22 and 60 T, and define B^ 1 as the av- 
erage inverse field taken over the width of a window. For 
the low field data set the Fourier window was 12.5-16 T. 

For both the a and /3 frequencies good quality fits are 
obtained for high and low angles, with effective masses 



2.35 ± 0.05 and m* = 1.22 ± 0.05 at 



10° 



expressed in units of m e and with error bars set by un- 
certainties in the least squares fit. The top panel shows 
data taken in the superconducting magnet system, giving 
a value for m* of 2.0 ± 0.3 at 9 = 0°, which is in agree- 
ment with the high field data, serving as a useful check 
of consistency between measurement techniques and sys- 
tems. In panel (c) data at higher angles is shown, demon- 
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FIG. 3: (color online) (a) The upper panel shows frequency 
(F) and effective mass (m*) as a function of angle for the 
a electron pocket. The lines are fits assuming conventional 
orbits about an ellipsoid formed from a parabolic band. The 
lower panel shows the same data for the 13 orbit, (b) Cartoon 
showing extremal orbits of the Fermi surface near the Dirac 
point arising from Eq. 3l k y , and k x lie in the plane of the 
orbit, while k z lies along the c axis of the reciprocal lattice. 
The j3 (7) orbit originates from the larger (smaller) pocket, 
while the a orbit encompasses a separate Fermi pocket derived 
from parabolic bands (not shown). 



strating an increase in m* for both orbits [IE] . 

We collect m* and F for all angles together in Fig. 3, 
with fitting errors reflecting the signal to noise, which is 
limited at higher angles by the shorter field range of the 
oscillations. The top panel shows data for the a orbit, 
which shows that the increase in F a is closely mirrored 
by a concomitant increase in m* . Intriguingly, the same 
cannot be said for the /3 orbit in the bottom panel, where 
Fg rises faster with angle than Trig. It is from these obser- 
vations that the presence of Dirac fermions in SrFe2As2 
can be inferred. 

To interpret our data we now consider two different 
scenarios. First, that the small Fermi surface pockets 
are ellipsoids derived from conventional parabolic bands, 
and second, that the pockets originate from bands with 
a linear dispersion near Ep. In the first case, we would 
expect the frequency as a function of angle to be given 
by BH: 



F(6) = 



m fi 



eh\J cos 2 9 + t; 2 sin 2 9 



(1) 



Where \x is the chemical potential, £ is the ratio of 
major to minor axes of the extremal elliptical cross sec- 
tion, and ?71q is the effective mass measured at 9 = 0. 



m*(9) has the same 9 dependence as F(9), with m*{9) = 
(eh/fj,)F(9). Fig. 3a shows the results of a best fit to this 
model. While a fit with parameters with /1 = 21 meV, 
itiq=2.2 and £=0.58 captures the angular dependence of 
the frequency and mass for the a orbit, it fails to explain 
the experimental data for the j3 orbit. The best fit to 
Fg(6) yields \x = 15 meV, m*=1.2 and £=0.25, which 
rises too fast with angle to account for the angular de- 
pendence of m* B . We are thus forced to look beyond a 
conventional model to describe this pocket. 

In the 122-pnictide materials, Dirac points are pro- 
posed to arise due to degenerate electron and hole bands 
that meet at Ep . The existence of a small amount of in- 
terlayer hopping provides access to these points, by cre- 
ating a series of stacked pockets along the fc^-axis that 
touch at their extrema. The influence of massless charge 
carriers originating from these regions of the Fermi sur- 
face on quantum oscillatory phenomena has been cap- 
tured in a minimal model |12j . which predicts a dip in 
the ratio m*{9)/F{9) at higher angles consistent with 
our observation for the /3 orbit. The starting point is a 
dispersion relation of the form: 



e = ±Hv*\k\ + 2tcos (ck z /2) + fj, 



(2) 



where c is the bilayer spacing for the body centered 
tetragonal crystal structure (12.3 A), v* is the charac- 
teristic Dirac velocity, t is the interlayer hopping pa- 
rameter, |fc| 



lk% + ky, and \x is the chemical poten- 
tial. For simplicity, the characteristic velocity v* is taken 
to be isotropic, although there is evidence for a small 
anisotropy from angle-resolved photoemission (ARPES) 
measurements [H]. As the angle 9 between the c-axis 
and the applied magnetic field is varied, the extremal 
orbits take approximately elliptical forms. Defining k x 
and k' as vectors lying in the planes of these orbits, as 
illustrated in Fig. 3b, gives 
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8 (3) 



This allows calculation of the expected angular vari- 
ation of the cyclotron effective mass and extremal cross 
sectional areas A^.e, which are directly related to the 
quantum oscillation frequencies via the Onsager relation, 
F e = (H/2-Ke)A ki e. 

Taking the ratio m*{9)/F{9) provides a straightfor- 
ward test of the predictions of this model. For Fermi 
surface pockets originating from bands with parabolic 
dispersion, we would expect a ratio which does not vary 
with angle, as has for example been shown in 2D organic 
metals [5D]. For bands with the dispersion parameter- 
ized by Eq. [2] a striking effect occurs - as 9 is increased 
the orbits begin to encompass regions of the Fermi sur- 
face close to the Dirac point, where velocities are high 
and masses very light [T7]- Thus, the orbitally averaged 
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FIG. 4: (color online) Ratio of the effective mass to cyclotron 
frequency as a function of angle for the a and /3 orbits. The 
solid line shows the fit to a model arising from the Dirac-like 
dispersion relation in Eq. [2] while the dotted line shows the 
expected behavior for conventional ellipsoids. 



cyclotron masses are smaller than those expected from 
parabolic bands, and the ratio m* / F is correspondingly 
reduced, providing an explanation for observation of the 
same. 

Figure 4 confirms this effect for the /3 orbit. In panel 
(a) m* a /F a is shown to be nearly constant over a wide 
range of angles, within the errors of the experiment, as 
expected for a Fermi surface with conventional parabolic 



bands. Panel (b) shows the same ratio for the f3 orbit, 
where a dip in m%/Fp is clearly observed, with the solid 
line representing a fit to a Fermi surface parameterized 
by Eq(3]with w* = 6.75 x 10 4 ins" 1 ,* = 15.2 meV and 
fi«0 meV. We conclude that this is direct proof of the 
presence of Dirac nodes sampled by the f3 orbit [21] . 

Despite the relative simplicity of the Fermi surface 
model derived from Eq(3| the extracted vt agrees rather 
well with that measurea with ARPES [THJ H2] , where a 
linear dispersion of bands in the vicinity of the Brillouin 
zone M point is observed. Through our observation of 
the angular dependence of quantum oscillations, we can 
unambiguously confirm the observation of linearly dis- 
persing bands that cross the Fermi energy near a Dirac 
node. Recent evidence from magnetotransport measure- 
ments on the same material interpret a linear magnc- 
toresistance at intermediate fields as evidence of Dirac 
excitations [33], which is also consistent with our find- 
ings. 

Having established the presence of Dirac nodes in the 
parent compounds of pnictide high T c superconductors, 
an interesting question remains as to the role they play 
in superconductivity in this system, and their influence 
across the doping or pressure phase diagram. Recent the- 
oretical work has suggested that the physical symmetry 
of the 122-pnictides and the topology of their band struc- 
ture lead to an unusual robustness of the nodal spin den- 
sity wave state [13] , with important consequences for the 
nature of the unconventional superconductivity induced 
by pressure or doping. 
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